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A generalized constant coupling approximation for quantum geometrically frustrated an- 
tiferromagnets is presented. Starting from a frustrated unit, we introduce the interactions 
with the surrounding units in terms of an internal effective field which is fixed by a self 
consistency condition. Results for the static magnetic susceptibility and specific heat are 
compared with previous results in the framework of this same model for the classical limit. 
The range of applicability of the model is discussed. 

PACS numbers: 75.10.Jm, 75.30.Cr, 75.40.Cx 

I. INTRODUCTION 

In the last several years, geometrically frustrated antiferromagnets (GFAF) have emerged as a 
new class of magnetic materials with uncommon physical properties, and have received a great 
deal of attention (see Refs. [T],^^ and references therein). In these materials, the elementary unit of 
the magnetic structure is the triangle, which makes it impossible to satisfy all the antiferromagnetic 
bonds at the same time, with the result of a macroscopically degenerate ground state. Examples 
of GFAF are the pyrochlore and the kagome lattices. In the former, the magnetic ions occupy the 
corners of a 3D arrangement of corner sharing tetrahedra; in the later, the magnetic ions occupy 
the corners of a 2D arrangement of corner sharing triangles (see Fig. p. In the case of materials 
which crystallize in the pyrochlore structure, the static magnetic susceptibility follows the Curie- 
Weiss law down to temperatures well below the Curie-Weiss temperature. At this point, usually 
one to two orders of magnitude smaller than the Need point predicted by the standard mean field 
(MF) theory, some systems exhibit some kind of long range order (LRO), whereas others show a 
transition to a spin glass state (SG). This is a striking feature for a system with only a marginal 
amount of disorder. Finally, there are some pyrochlores which do not exhibit any form of order 
whatsoever, and are usually regarded as spin liquids. In the case of the kagome lattice, even though 
there are very few real systems where this structure is realized, the magnetic properties fall in two 
major categories: the vast majority of the compounds studied show a transition to a LRO state with 
a non collinear configuration of spins, and a few systems exhibit no LRO, but a SG like transition. 

For these reasons, it is easy to understand the large amount of attention, both from the experi- 
mental and the theoretical points of view, these systems have attracted. From the theoretical point 
of view, a number of techniques have been used to try to understand the origin of the puzzling 
properties mentioned above. All the theoretical results seem to indicate that the classical Heisen- 
berg model with only nearest neighbor interactions does not dispkv anvlong range order for these 
geometries, in agreement with Monte Carlo (MC) result ^fflfflBBElli. 

There are also a relatively few works which have dealt with the quantum effects in these 
systemsolBnOliSoEltilEl, even though the main interest during last years has been on the classi- 
cal GFAF. However, in a recent worku, the present authors showed that there are some features of 
the experimental data for the susceptibility in pyrochlore compounds that can only be understood 



(a) Pyrochlore lattice. 



(b) Kagome lattice. 



FIG. 1: The magnetic lattices considered in this work. 

in a quantum framework, as, for example, the maxima appearing in this quantity at temperatures 
well below the Curie-Weiss temperature. 

In another recent workS, the present authors developed a generalization of the well known 
constant coupling method0 (CC), which can be applied to frustrated geometries, the so called 
generalized constant coupling (GCC) method. This technique was applied successfully to the 
study of the classical Heisenberg Hamiltonian with nearest neighbors (NN) interactions in both 
the pyrochlore and kagome lattices. In spite of the mathematical simplicity of this technique, the 
results obtained for the susceptibility are essentially exact when compared with Monte Carlo data 
in both types of lattices. Moreover, the calculated specific heat is also in very good agreement with 
Monte Carlo data, even though there are some deviations at very small temperatures, which can 
be understood in the light of the fact that a MF theory cannot properly describe this property (or, 
equivalently, the internal energy) at T = K, due to the distinct nature of the excitations which are 
important in this limit, namely, the spin waves. In any case, the GCC method provides an excellent 
starting point for studying the magnetic properties of the frustrated systems in the paramagnetic 
region. 

The first question that arises is which features of these magnetic properties are different in 
a quantum formulation of the GCC method, and that is precisely what we try to answer in this 
work. However, in contrast with the classical limit, there are no reliable quantum Monte Carlo 
calculations for these systems, due to well known difficulties that arise in the quantum version of 
this technique, which makes it very difficult to check the accuracy of the present quantum version 
of the GCC method. It is tempting to compare the results of this method with experimental data 
available in the literature. However, here we are considering the most restrictive approach to the 
problem by focusing on Heisenberg models with NN interactions only, whereas in real systems, 
as pointed out above, there are always additional effects, so we think such a comparison makes no 
sense at this point. Therefore, the most we can do is to give the results for the quantum version 
of the GCC method, and study how they evolve towards the classical results for large values of 
the individual spin quantum number. If this transition from the quantum to the classical limit is 
smooth, we can expect the predictions of the model to be, at least, reasonable in the paramagnetic 
region. 
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TABLE I: Values of g(S) for some representative values of s . Integer values of S correspond to 
the pyrochlore case (p — 4) for any value of s and the kagome lattice for integer values of s . 
Half integer values of S correspond to the kagome lattice (p = 3) for half integer values of s . 



H. THE MODEL 

The Heisenberg model with only nearest neighbor INN) interactions in the presence of a mag- 
netic applied field H is described by the HamiltonianH 

H = J^Si ■ Sj - H ^2s Zi , (1) 

(m) * 

where J is the positive antiferromagnetic coupling, s, and Sj represent quantum spins of modulus 
s located in a pyrochlore or kagome lattice and s Zi the corresponding component along the applied 
field, and the sum is done over pairs of NN. 

The idea of our approximate method is based on the experimental fact that the spin-spin corre- 
lations in the GFAF lattices are short rangedP. Therefore, it is a reasonable approximation to start 
by considering isolated units (tetrahedra or triangles, for the pyrochlore and kagome lattices, re- 
spectively), and later add the interactions with the surrounding units in an approximate way. Thus, 
it is important to first study the properties of the individual units. This task has been carried out 
by Garcfa-Adeva and Huber for the quantum casaJj, and we will not repeat the derivation here. 
It will suffice to say that the susceptibility per spin of an isolated unit with p spins (p = 4 for the 
pyrochlore and p = 3 for the kagome) is given by 



v m = vg> _ 1 J2s9(S)S(S+l)(2S+l)e-^^ 
Xp{ ' 3pT 3pT J2s9( S )( 2S + 1 ) e ~ jS{S+1)/2 



where j = J/T, and (S 2 p ) represents the average value of the square of the total spin of the unit, 
and g(S) is the number of configurations with total spin S (see Table |). 

For some derivations, it is better to express all the thermodynamic quantities in terms of the 
function 



(s 2 v ) 

e P (T) = P ' u - 1. (3) 
ps {s + 1) 
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Therefore, the susceptibility of the isolated unit can be expressed as 



XriT) = *(* + i)(i + *,(*■» (4) 



It is important to notice that the only difference from the derivation carried out in Ref. is in 
definition of (Si*). The rest of the method is essentially the same. 

Next, the interaction with neighboring units is introduced as an unknown internal effective field, 
Hi, created by the (p — 1) NN ions outside the unit. The CC approximation consists of taking this 
internal field as Hi = (p — l)H', where H' is the average internal field acting on a spin due to 
each of its NN. By comparison, in the case of the standard CW model, the internal effective field 
is given by Hi = 2(p — 1)H', as the ions are considered separately, and each has 2(p — 1) NN 
in the corner sharing structures considered in this work (see fig. p. The internal field is evaluated 
by imposing the self consistent condition of equating the magnetization per spin in the field with 
that of a unit in the field. In the paramagnetic region, and small applied field, the susceptibility per 
spin is then given byS 

v gcc (T] = sp(s + 1) 1 + £ p (T) 

Xp y ' 3T l-e p (T)' K J 

The corresponding expressions for the internal energy0and the specific heat per spin are given by 

u p (T) = Js (s + l)eJT) (6) 



and 



respectively. 



c p (T) = J So ( So + l)A £p ( T ) ; ( 7) 



III. RESULTS FOR THE SUSCEPTIBILITY AND THE SPECIFIC HEAT 

In this section, we will present the main results of the model introduced in the previous sec- 
tion. One is tempted to compare them with the available experimental data. However, such a 
comparison makes no sense at this point because, as stated in the introduction, there are several 
effects not included in this simplified model such as dipole-dipole interactions, anisotropics, or 
further neighbor interactions, to cite some examples. Therefore, we will present the results of the 
present model and study how the quantum case evolves towards the classical limit for relatively 
large values of the individual spins. If the quantum behavior of the physical quantities of interest 
(namely, susceptibility and specific heat) is not so different from the classical one for relatively 
large values of sq, we can feel confident that the predictions are essentially correct, at least in the 
temperature region where deviations from the classical behavior are small, as we know that the 
GCC model gives an accurate description of these quantities in whole the temperature range. It 
would be desirable also to compare the predictions of this model with the ones obtained from more 
sophisticated techniques (quantum Monte Carlo, high temperature series expansions, density ma- 
trix renormalization group, exact diagonalization, and son on). However, that comparison is not 
possible in general, due to the geometrical complexity of GFAF which prevents one from applying 
any of these techniques to arbitrary values of s , and these kind of calculations have been only 
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(a) Magnetic susceptibility. (b) Specific heat. 



FIG. 2: Susceptibility and specific heat for the pyrochlore lattice for various values of s : (long 
dashed line ) so — 1/2; (short dashed line) s = 1; (dotted line) so = 3/2; (dotted dashed line) 
so = 7/2 . The solid line is the classical limit of the model. The double dotted line in the case of 
the susceptibility represents the Curie-Weiss law. 

carried out for the quantum s = 1/2 which, as we will see below, is the case which we expect to 
have the worst results. 

Let us start by considering the results for the susceptibility and specific heat in the pyrochlore 
lattice, which can be seen in Fig. £2] for some values of s . As we can see from the observation of 
that figure, the classical susceptibility deviates from the Curie-Weiss behavior below Qcw- The 
quantum case is qualitatively similar to the classical behavior, except at very low temperatures, 
where it passes through a maximum and then falls to zero. However, we can also see that even 
for relatively small values of sq the quantum behavior is not so different from the classical one 
down to the maximum. Regarding the predictions for the specific heat, again, the main difference 
is that it goes through a maximum, and later falls to zero, in contrast with the classical situation, 
where it goes to a constant value at T = K. Additional information can be gained by studying 
the dependence of the position of the maxima with the value of s m Fig. |3|. In that figure can be 
seen how nicely both the maxima in the susceptibility and specific heat go to the zero classical 
value. 

Let us now turn to the kagome lattice. In this case, it is reasonable to expect the model to 
work less well than in the pyrochlore lattice, as it is well known that the role of long wavelength 
fluctuations is more important in 2D lattices, which are explicitly neglected in any MF approach. 
However, from the MC results for the classical limit, it turns out that long wavelength thermal 
fluctuations are unimportant for GFAF at finite temperatures, which partially explains why the 
classical limit of our model is essentially exact when compared to MC data for the susceptibility. 
The quantum case is a more delicate one, as quantum fluctuations can play an important role. 
However, these fluctuations will be only important at low temperatures, in a region where our 
model does not apply for other reasons we will explain below. 
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(a) Pyrochlore lattice. (•) Maxima in the sus- 
ceptibility; (o) Maxima in the specific heat. 



FIG. 3: Position of the maxima appearing in the 
guide to the eye. 
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(b) Kagome lattice. (□) Maxima in the 
susceptibility; (o) Maxima in the specific 
heat for half integer so; (x) Maxima in 
the specific heat for integer values of sq. 

sceptibility and specific heat. Solid lines are a 



In order to proceed further, we will split the discussion for the kagome lattice in two cases, 
corresponding to integer and half integer values of the individual spins. 

Let us first analyze the results for integer values of s > which can be seen in Fig. |j. There are no 
remarkable differences with respect to the discussion for the pyrochlore, except for the fact that the 
maximum in the susceptibility is present even in the classical limit. The only qualitative difference 
is the later fall to zero of both the susceptibility and specific heat. Also, we have computed the 
evolution of the maxima for both the susceptibility and specific heat. Again, both quantities go 
smoothly to their classical values. It is important to stress that the fact that the susceptibility goes 
to zero is related to the ground state total spin of the cluster, which, for integer values of the 
individual spins, is always a singlet. This is not the case for half integer values of s . 

For half integer values of s , there are very important differences with respect to the previous 
cases, as is evident in Fig. The susceptibility goes through a maximum, then reaches a minimum, 
and diverges as T — > 0. This is a consequence of the fact that the ground state of the triangular 
cluster for half integer values of s is a doublet in our model. However, it is now generally accepted 
that the real ground state of the kagome lattice is a singlet, and may have a gap in the spectrum so, 
actually, the susceptibility should either be finite or go to zero as T — > 0. Obviously, the present 
breakdown at very low temperatures is not a particular feature of the mean field approach, but it is 
shared by every model which starts from finite cluster calculations, with an odd number of spins 
in the elementary cluster. Moreover, the maxima in the specific heat follow a different trend than 
in the integer s case, even though, the classical limit is correct (see Fig. 0). Even more striking 
is the fact that both the susceptibility and specific heat for so = 1/2 are qualitatively different 
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(a) Magnetic susceptibility. (b) Specific heat. 



FIG. 4: Susceptibility and specific heat for the kagome lattice with integer values of s : (long 
dashed line) s = 1; (short dashed line) s Q = 2; (dotted line) s = 3; (dotted dashed line) s = 4. 
The solid line is the classical limit of the model. The double dotted line in the case of the 
susceptibility represents the Curie-Weiss law. 



from all the previous presented results. The susceptibility does not show a maximum, but diverges 
as 1/T as T approaches 0. The maximum in the specific heat does not follow the extrapolated 
behavior for other half integer values of sq. These kinds of problems have been already pointed 
out in the literature for more rigorous methods than ours. A qualitative argument that can help 
to understand these deviations is the fact that s = 1/2 is the "most quantum" case, in which the 
physical quantities are more sensitive to the discrete structure of the energy levels of the system. 
Taking in to account that for half integer values of s , the energy level structure of our model is 
qualitatively incorrect, we should not be surprised by this special behavior. 

As commented above, if we want to further check the range of applicability of the model, it 
would be desirable to compare its predictions with the ones obtained from more sophisticated 
methods, as exact diagonalization of small clusters or high temperature series expansions. Unfor- 
tunately, these calculations have been carried out only for the kagome with sq = 1/ 200 where, as 
we have argued above, we can expect the model to be least successful. However, even though we 
do not expect complete agreement, we have carried out such a comparison with results obtained 
from the aforementioned methods and the results are presented in Fig. ^. In that figure we have 
also plotted the corresponding curves corresponding to s = 1 and s = 3/2 for comparison, at a 
qualitative level. It is important to stress that the curves depicted there do not contain any fitting 
parameter nor have they been rescaled in any sense. Surprisingly, the susceptibility calculated 
from our model for so = 1/2 is in very good agreement down to ~ 0.5 Qcw with both high tem- 
perature series expansions and the diagonalization of a cluster with 18 spins. The corresponding 
curve for the specific heat is in poor quantitative agreement with the results from high temperature 
series expansions. However, it describes accurately the position of the maximum. In any case, 
we can expect the agreement with the specific heat to be worse, as it considers independent units, 
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(a) Susceptibility per ion. 
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FIG. 5: Susceptibility and specific heat for the kagome lattice with half integer values of s : (long 
dashed line) s = 1/2; (short dashed line) s = 3/2; (dotted line) s = 5/2; (dotted dashed line) 
so = 7/2 . The solid line is the classical limit of the model. The double dotted line in the case of 
the susceptibility represents the Curie-Weiss law. 



neglecting the emerging correlations at low temperatures. Moreover, we can see that the calcu- 
lated susceptibilities for So = 1 and s = 3/2 reproduce the main qualitative features of the high 
temperature expansion and finite cluster results. At this point, it is difficult to say if this a fortunate 
coincidence or an indication of the accuracy of our model for higher values of s . Numerical data 
for these values of the spins would be necessary to decide between the two cases. 

In any case, it is important to notice that some of the qualitative features obtained in this sim- 
ple model, have been observed in experimental studies. Again, it is difficult to say at this point 
if those features are due to additional interactions not accounted for in this model. For example, 
the existence of a maximum in the susceptibility for the pyrochlore lattice has been systemati- 
cally observed in all the experiments at low enough temperatures. Actually, these maxima where 
successfully interpreted by the present authors with an even simpler model, in which the interac- 
tions with nearest neighbors and next nearest neighbors were taken into accountEI. Moreover, in 
the work by Wills and coworkers^ on iron jarosites, which is one of the few systems where the 
kagome lattice has an experimental realization, the magnetic susceptibility exhibits the maximum 
and a later upturn, at a temperature comparable with that predicted by our model. However, in 
that temperature region, it seems likely that the dilution of the magnetic lattice by non magnetic 
impurities plays an important role, giving rise to some kind of spin glass behavior, for which it is 
very difficult to extract any conclusion. 

In the light of these results, we think that a prudent estimate of the breakdown of our model is 
set by the position of the maximum in the specific heat. As we approach this temperature, emerg- 
ing correlations that cannot be described in the GCC formalism enter into play. For the pyrochlore 
lattice and the kagome with integer values of s , it is very likely that the present description is 
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FIG. 6: Comparison of the GCC predictions with results from finite size clusters calculations and 
high temperature series for the kagome 1/2 lattice^: (solid line) GCC model for s = 1/2; (long 
dashed line) GCC model for s = 1; (dashed line) GCC model for s = 3/2; (o) Results from 
high temperature series expansions; (•) Results from exact diagonalization of a cluster with 18 
spins. The double dotted line in the case of the susceptibility represents the Curie-Weiss law. 

qualitatively correct for temperatures well below the peak temperature of the specific heat. How- 
ever, we think that the description below that temperature for half integer s in the kagome lattice 
is incorrect. Probably, there is a maximum and a minimum in the susceptibility, but we expect 
the susceptibility either vanishes or goes to a finite value at T — > 0, not to diverge. In any case, 
it is still a significant improvement with respect to the Curie- Weiss theory and other MF theories. 
However, more rigorous techniques are needed in order to verify these assertions. 



IV. CONCLUSIONS 

In this work we have presented the quantum version of the generalized constant coupling model, 
which was shown by the authors to be in excellent agreement with Monte Carlo data for the 
susceptibility and specific heat in the classical limit, for both the pyrochlore and kagome lattices. 
There are some important qualitative differences between the classical and quantum behaviors 
which are important at low temperatures. 

The main features of the susceptibility in this model are the presence of maxima in both the 
susceptibility and specific heat, similar to those found experimentally. The susceptibility for the 
kagome lattice with half integer values of the individual spins is found to pass through a maximum 
and then, after a minimum, diverges as T — > 0. This divergence is due to a non zero value of the 
total spin of the ground state of the elementary units for these values of the individual spins. There 
is some experimental evidence for the existence of such an upturn in the susceptibility in the iron 
jarosite systems. However, the divergence of this quantity in our model, which is due to the finite 
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size cluster with odd number of spins, is incorrect. 

The results for the kagome 1/2 are compared with high temperature series expansions and exact 
diagonalization of small clusters. The susceptibility is found to be in good agreement with those 
results down to T ~ 0.5 ®cw, which is a remarkable fact for so simple a model. Even though 
the specific heat is not in quantitative agreement at this temperature, the position of the maximum 
is adequately predicted. Results from the aforementioned techniques for these lattices for higher 
values of so would be desirable in order to check the accuracy of the GCC model. 

In any case, we feel that the present model provides an adequate description of the nearest- 
neighbor kagome and pyrochlore systems down to the temperature of the peak in the specific heat. 
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